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Abstract: A new hybrid experiment has been started by AS7 experiment at Tibet, China, since August 2011, which 
consists of a low threshold burst-detector-grid (YAC-11, Yangbajing Air shower Core array), the Tibet air-shower array 
(Tibet-Ill) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement 
of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-11) is investigated by means of an 
extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-11) array 
and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the 
new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of 
the major component at the knee. 
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1 Introduction 

The all-particle energy spectrum of primary cosmic rays 
is well discriminated by a power law dN/dE ex E -7 over 
many orders of magnitude, with 7 changes sharply from 
2.7 to 3.1 at about 4 PeV [1]. The break of the all-particle 
energy spectrum is called the "knee", and the correspond- 
ing energy range is called "knee region". Although the 
existence of the knee is well confirmed by experiments, 
there still is controversy on its origin. In order to explain 
the sharpness of the knee, two scenarios [2] (model A and 
model B) are proposed as shown in Fig.l and Fig. 2. 

In model A, an excess component is assumed to overlap 
the global component, and its spectrum shape suggests that 
it can be attributed to nearby source(s) because it is sur- 
prisingly close to the expected source spectrum of the dif- 
fuse shock acceleration. Middle composition is predicted 
by this model at the knee. In model B, a hard observed 
energy spectrum of each element from a given source is 
assumed. The sharp knee can be explained by a rigidity- 
dependent acceleration limit and the hard spectrum due to 
nonlinear effects. Iron-dominant composition is predicted 
by this model at the knee and beyond. In order to distin- 
guish between Model A and Model B and many other mod- 
els, measurements of the chemical composition around the 
knee, especially measurements of the primary spectra of 
individual component till their knee will be essentially im- 
portant. Therefore, we have developed the Yangbajing Air 
shower Core detector (YAC) which is operated along with 
Tibet air- shower array (Tibet-Ill) and underground water 
Cherenkov muon detector array (MD) simultaneously, as 
shown in Fig. 3, and the second phase of YAC is so called 
YAC-II. In this paper, the capability of the measurement of 
the chemical components with use of the (Tibet-III+ YAC- 
II) is investigated. The simulation results using (Tibet- 
III+YAC-II+MD) will be reported in the near further. 

2 A New Tibet Hybrid Experiment 

The Tibet new hybrid experiment (Tibet-III+ YAC-II) has 
been operated in Tibet, China, since August, 1st, 201 1. The 
merit of this experiment is that the atmosphere depth of the 
experimental site (4300 m above sea level; 606 g/cm 2 ) is 
close to the maximum of the air shower development with 
energies around the knee and the shower maximum val- 
ues are almost independent of the masses of primary cos- 
mic rays. The Tibet-Ill consists of 789 detectors while the 
YAC-II consists of 124 detectors (as shown in Fig. 3). The 
inner 100 plastic scintillator units of YAC-II are arranged as 
an array (10x10 grid) each with an area of 50 cm x 80 cm, 
with 1.875 m interval; and the outer 24 plastic scintillator 
units are arranged around the inner array each with an area 
of 50 cm x 100 cm. The outer 24 units are used to reject 
non core events whose shower cores are far from the YAC- 
II array. Each detector of YAC-II consists of lead plates 
with a thickness of 3.5 cm above the scintillator to con- 
vert high energy electrons and gammas to electromagnetic 
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Figure 1 : Model A: Sharp knee is attributed to extra com- 
ponent from nearby source. 
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Figure 2: Model B: Nonlinear effect in the diffusive shock 
acceleration mechanism can explain the structure of the 
knee when high acceleration efficiency for heavy elements 
is assumed. 
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Figure 3: The Tibet-III+YAC-II+MD array. The Tibet-Ill 
consists of 789 detector units, the YAC-II consists of 124 
detector units. The five MDs in the blue frame are set up 
this year and acquiring data soon. 



showers. Each unit of YAC-II is attached with two photo- 
multipliers of high-gain (HAMAMATSU: R4125) and low- 
gain (HAMAMATSU: R5325) to cover the wide dynamic 
range from 1 MIP (Minimum Ionization Particle) to 10 6 
MIPs. The hardware of YAC-II is described in [3]. 

The area of YAC-II array is about 500 m 2 , locating 
near the center of the Tibet-Ill, and operating simultane- 
ously with it. For an air shower event, the Tibet-Ill pro- 
vides the arrival direction (0), the air shower age (Age) and 
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the air shower size (N e ) which is interrelated to primary en- 
ergy, the YAC-II measures the high energy electromagnetic 
particles in the core region so as to obtain the characteris- 
tic parameters of air- shower cores. When a YAC event is 
triggered, its accompanying air shower is simultaneously 
recorded. The matching between an AS and a YAC event 
is made by their arrival time stamps. The air-shower direc- 
tion can be estimated with an error smaller than 0.1° above 
100 TeV, and the primary energy resolution is estimated to 
be 12% at energies around 10 15 eV by our simulation. 
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3 Simulation and Analysis 

A Monte Carlo (MC) simulation has been carried out 
on the development of extensive air showers in the atmo- 
sphere and the response in YAC-II and Tibet-Ill. The simu- 
lation code CORSIKA (version 6.024) including QGSJET2 
hadronic interaction models [4] is used to generate AS 
events. Considering the aim of this simulation is just to 
check the capability of the hybrid experiment, QGSJET2 
interaction model and heavy dominated (HD) [1] primary 
composition model are used. Primaries isotropically inci- 
dent at the top of the atmosphere within the zenith angles 
from to 60 degrees are injected into the atmosphere. The 
simulated data were analyzed in the same manner as in the 
procedure for the experimental data analysis. The electro- 
magnetic showers in the lead layer induced by electrons or 
photons that hit any detector unit of the array are simulated 
by a subroutine which is based on the detector simulation 
code EPICS [5]. The detector performance, the trigger ef- 
ficiency of detectors and the effective area are adequately 
taken into account based on the experimental conditions. 
Normally, the following quantities of YAC-II are used to 
characterize an air- shower core event: 

N b - the number of shower particles under the lead plate of 
a detector unit; 

Nhit - the number of "fired" detector units with N b > a. 
given threshold value; 

N b top - the maximum burst size among fired detectors; 

^N b - the total burst size of all fired detector units; 

< R > - the mean lateral spread, < R >=%2ri/(Nhit-l)', 

<N b R> - the mean energy-flow spread, 

<Nj ) R>=^2(Nj ) iXri)/Nhit, where N bi and n are the burst 
size in the i th fired detector unit and the lateral distance 
from the air shower core to the center of the i th fired detec- 
tor, respectively. 

We divided the MC data into two datasets. Due to the dif- 
ference of the detection efficiency, the first dataset selects 
events that are enriched with proton+helium origin (called 
tagged-I dataset), while the second dataset contains events 
of heavy-primary origin (called tagged-II dataset). The fi- 
nal select conditions for tagged-I and tagged-II are as fol- 
lows: 

(1) N b > 100, N hit > 12, N h top > 3000, N e > 1 x 10 5 ; 

(2) N b > 100, N hit > 20, N h top > 1500, N e > 5x 10 5 . 



Figure 4: The detection efficiency SQA e ff of proton, pro- 
ton+helium and iron. 

Besides, the detector unit with N b top is requested to be 
located at inner 8x8 grid for both data-sets. Then, we ob- 
tain 1.43 xlO 5 and 5.17xl0 4 events for the tagged-I and 
tagged-II dataset, respectively. The detection efficiency 
SQA e ff is shown in Fig.4 . 

The separation of the primary mass is realized by use of 
a feed-forward artificial neural network (ANN[6]) method, 
whose applicability to our experiment was well confirmed 
by the MC simulation [7]. For one thing, we need to sepa- 
rate protons from other nuclei by training the network with 
a proton flag, and then separate proton+helium from other 
nuclei by training the network with a proton+helium flag, 
thus we can get the helium energy spectrum by subtract- 
ing the number of proton events from the proton+helium 
events, so does iron. The following 8 parameters are 
input to the ANN with 40 hidden nodes and 1 output 
unit: (l)N hi u (2)N b top , (3)£JV 6 , (4)< R >, (5)<N b R>, 
(6)N e , (7)0, (8) Age, where the first five parameters are 
given by YAC-II, and the last three are obtained by Tibet- 
Ill. 

4 Results and Discussion 

The ANN training results of proton, proton+helium and 
iron are presented in Fig. 5, where average purity and se- 
lection efficiency over whole energy range are shown. The 
ANN output value T is used to separate the primary nuclei 
groups. In this paper, the events with T < T c (or T > T c 
) are regarded as proton or proton+helium group (or iron 
group). The value T c , purity (p) and selection rate (e) of 
events satisfying the criterion at various energy regions are 
summarized in Table 1 . 

The Fig. 6 shows the estimated primary energy spectra of 
proton, helium and iron compared with the assumed ones. 
One can see that the assumed primary energy spectrum of 
proton, helium and iron are well reproduced by the esti- 
mated ones respectively, and the estimated ones could well 
connect with the results obtained by direct observation, as 
shown in Fig. 7. It needs to be remarked that the iron spec- 
tra seem to be higher than the observed ones, just because 
of the HD model we used. The results above show that the 
new burst detector array is powerful to study the chemi- 
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Table 1: The purity (p)(%) and the selection rate (e)(%) of the selected primary groups. 
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Figure 6: Energy spectra of primary proton, helium and 
iron. 
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Figure 7: The expected spectra of proton, helium and iron 
obtained by MC compared with other experimental data. 



Figure 5: The ANN training results of proton, pro- 
ton+helium and iron. The average purity and selection rate 
over whole energy range of protons are 88%, 46% at T c = 
0.2, while 96%, 28% for proton+helium at T c = 0.1, 80%, 
54% for iron at T c = 0.7. 

cal composition, in particular, to obtain the primary energy 
spectrum of the major component at the knee. 
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